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Abstract 

First principles total-energy pseudopotential calculations have been performed to investigate 
STM images of the (110) cross-sectional surface of Mn-doped GaAs. We have considered configu- 
rations with Mn in interstitial positions in the uppermost surface layers with Mn surrounded by As 
(Int/is) or Ga (Intg-a) atoms. The introduction of Mn on the GaAs(llO) surface results in a strong 
local distortions in the underlying crystal lattice, with variations of interatomic distances up to 3% 
with respect to unrelaxed ones. In both cases, the surface electronic structure is half-metallic (or 
nearly half metallic) and it strongly depends on the local Mn environment. The nearby Mn atoms 
show an induced spin-polarization resulting in a ferromagnetic Mn-As and antiferromagnetic Mn- 
Ga configuration. The simulation of the STM images show very different pattern of the imaged 
Mn atom, suggesting that they could be easily discerned by STM analysis. 

PACS numbers: PACS: 73.20.At; 75.50.Pp; 75.70.Rf; 71.55.Eq 
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I. INTRODUCTION 



The easy integration of ferromagnetism with semiconducting properties in the same host 
material provided by Diluted Magnetic Semiconductors (DMS's) has been considered an 
important breakthrough in the semiconductor microelectronics. This is mainly due to the 
unprecedented opportunity to create a new class of device which would combine the spin 
degree of freedom to process, to transfer as well as to store information. Spintronics is the 
emergent technology which exploits the quantum propensity of electrons to spin as well as 
making use of their charge state.^^i^ 

The discovery of ferromagnetism in Mn-doped GaAs semiconductor has become a mile- 
stone in spintronic revolution: Mn^Gai_^As alloys are directly related to the existing GaAs 
technology, resulting in the practical realization of device structures combining ferromagnetic 
and nonmagnetic layers.— 

There are several possibilities for a single Mn to be incorporated in the GaAs. It can 
occupy or the cation site (substitutional Mn, Mn^a) either the anion site (As antisite, Mn^^); 
it can also occupy interstitial sites, as reported by K.M. Yu et al^ Further, other structural 
defects could be present in the alloy, such as As antisite (Asg^). The fraction of Mn dopants 
occupying one or another location depend on the growth conditions and techniques .— 

The Curie temperature (Tc) is a key parameter in designing room-temperature spintronic 
devices. The highest T^ reachable for Mn^-Gai-^^As up to few years ago was 110 K,— i.e. 
rather low for practical technological purposes. It has been shown that interstitial Mn atoms 
have a crucial role in magnetic properties of the samples.^^ An intense experimental and 
theoretical efforts have been pursued in the last years in order to understand the physics of 
this material and how to raise the Curie temperature. 

Nowadays, a new method has been proposed as alternative to the growth by Molecular 
Beam Epitaxy (MBE) of bulk Mn^-Gai^^As random alloy: the dopant atoms are incorpo- 
rated in the sample in such a way to give rise to a Dirac's 5 function concentration profile 
(with locally high dopant concentration) along the grow direction ((5-doping).— Remarkably, 
an important enhancement of Tc is obtained in these 5-doped samplesi^ (the highest Tc 
obtained so far with 5-doped sample is 250 K).^° Very recently, Mn J-doped GaAs samples 
in (001) direction have also been grown at TASC Laboratory in Trieste.— 

Therefore clarifying the site geometry and the local environment of impurities in 5-doped 



2 



GaAs:Mn should shed hght on the understanding and the optimization of the magnetic 
properties of the system. From the experimental point of view, this study can be pursued 
with cross-sectional Scanning Tunneling Microscopy (XSTM): the Mn-doped GaAs samples 
are cleaved along the natural (110) cleavage plane and then analyzed by STM microscopy. 

In recent years, several XSTM studies on Mua-Gai-^^As alloys have been performed but 
the local environment (and preferential geometric site) of defects has not been clarified 
yg^,i2ji3ji4ji5 pfQjji the theoretical point of view, the existing simulated XSTM images have 
mainly focused on the characterization of substitutional impurities on uppermost surface 
layers, while a complete and detailed investigation of interstitial impurity on uppermost 
surface layers is still lacking thus preventing the possibility of a full interpretation of the 
new XSTM images acquired. 

Therefore, stimulated by the recent growth and following XSTM analysis of Mn 5-doped 
GaAs samples at TASC,— we have performed density functional calculations to investigate 
the structural, electronic and magnetic properties of a single Mn dopant, by focusing our 
attention on the impurity interstitial surface configurations. We have also simulated the 
corresponding STM images. 

This paper is organized as follows: in the next section we describe the computational 
method; in Sect. 3 we present our results for the structural, electronic and magnetic prop- 
erties; in Sect. 4 we discuss our results for the XSTM images; finally, in Sect. 5 we draw 
our conclusions. 

II. COMPUTATIONAL DETAILS 

Our study has been performed within Density Functional Theory (DFT) framework in 
the Local Spin Density Approximation (LSDA) for the exchange-correlation (XC) functional 
by using state-of-the-art first-principles pseudopotential self-consistent calculations, as im- 
plemented in the ESPRESSO/PWscf code.^^ We used the scheme of Ceperley and Adlep^ 
(with the parametrization of Perdew and Zunger— ) for XC functional. Mn atom is de- 
scribed by an ultrasoft (US) pseudopotential (PP)-i^ while norm-conserving PPs have been 
considered for Ga, As and H atoms. 

Test calculations have shown that a kinetic energy cutoff for the wave functions equal to 
22 Ry and a 200 Ry cutoff for the charge density are sufficient to get well converged results. 
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We estimate the numerical uncertainty to be ~ 0.01 A for relative atomic displacements 
and ~ 0.02 //^ for the magnetic moments. The relaxed internal atomic positions have 
been obtained by total-energy and atomic-force minimization using the Hellmann-Feynman 
theorem.'^ 

The surface is modelled with periodically repeated cell containing one Mn atom; a (110) 
slab geometry with a 4x4 in-plane periodicity has been used. The simulation cells are made 
up of 5 atomic layers and a vacuum region equivalent to 8 atomic layers. The bottom layer 
has been passivated with Hydrogen atoms in order to simulate semi-infinite bulk material.— 
In the energy minimization only the three uppermost layers are allowed to relax, while the 
others are considered bulk-like. 

Two different configurations have been considered for Mn on the surface, namely lTitA_s{Ga) 
with As (Ga) atoms as nearest neighbor atoms. In each case, the distances between the Mn 
atom and its periodic image on the (110) plane are 15.7 A along the [110] and 22.2 A along 
[001]. 

XSTM images are obtained within Tersoff-Hamann model,— where the constant current 
STM images are simulated from electronic structure calculations by considering surfaces of 
constant integrated local density of states. 

III. STRUCTURAL, ELECTRONIC AND MAGNETIC PROPERTIES 
1. Structural properties 

The GaAs(llO) surface is well known from experimental as well as theoretical point of 
view.^^ In figO we show a ball and stick model of the clean surface, side and top views. 
The surface unit cell in shown in the top view. In this and the other figs., black spheres are 
cations (Ga atoms), grey spheres are anions (As atoms). 

At the top layer, the Ga surface atoms relax inward while the As atoms are shifted above 
the surface. Due to overbinding in the LDA approximation, our theoretical GaAs lattice 
constant (5.55 A) is smaller than the experimental one (5.65 A) but the relevant calculated 
structural parameters for the clean surface such as Ai^^ (relative displacement of the anion 
and cation positions in the uppermost layer, normal to the surface) and a (the buckling 
angle), shown in figO are 0.68 A and 30.36° respectively, which well compare with the 
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experimental values 0.65±0.03 A and 27 A°^^^ and other theoretical works.— i2Si21 

In zinc-blende bulk crystal there are two inequivalent tetrahedral interstitial positions 
for Mn which differ in their local environment: we denote them as Int^s or Int^a according 
whether Mn is surrounded by As or Ga atoms respectively. There is also an hexagonal 
interstitial position where Mn is surrounded by both As and Ga atoms. In figEl we show 
the different cases. The tetrahedral interstitial site in the ideal geometry has four nearest- 
neighbor (NN) atoms at a distance equal to the ideal host bond length di and six next- 
nearest-neighbor (NNN) atoms at the distance ^2 = which are Ga(As) atoms for 
IntAs(Ga)5 respectively. In the hexagonal interstitial position the Mn is surrounded by 3 
As and 3 Ga atoms at distance \f^di. Throughout this work we have considered only 
tetrahedral interstitial position (the total energy corresponding to the hexagonal interstitial 
site is higher by more than 0.5 eV) .■^'^^'^SiSi 

In fig. 121 we show a ball and stick side (a) and top (b) view of the relaxed Int^^, Intca 
configurations. Only the three topmost layers and the atoms closest to Mn are shown. 
Black spheres are cations (Ga atoms), grey spheres are anions (As atoms); Mn is explicitly 
indicated. In the relaxed structure, due to symmetry breaking because of the surface and 
the consequent buckling of the outermost surface layers, the NN and NNN bond lengths are 
no longer equal. Furthermore, some relaxed NNs bond lengths turn out to be longer than 
NNNs ones. In the following, we do not longer distinguish among NN and NNN (they are 
referred simply as NN atoms) but we simply refer to surface and subsurface atoms, as shown 
in the figure. 

The two relaxed configurations differ in energy by ~ 130 meV/Mn atom (Int^a is 
favoured). This is in contrast to the bulk case, where it has been found that they differ 
only by ~ 5 meV/Mn"^^ and Int^s instead is slightly favored. We have tested the reliability 
of our final relaxed interstitial configurations by considering different starting geometries 
(details in Ref. l32i ). other than the simple ideal (110) truncated bulk. In all cases, the final 
relaxed configuration is the same. 

The atoms with the most sizeable displacements from the ideal zinc blende positions are 
the Mn impurities and their neighbors, on surface or subsurface. In Tab. H] we report the 
inward/outward relaxations respect to the ideal (110) surface plane. 

In Int^s; Mn relaxes outward by ~ 0.06 A and Asgurf i^^subsurf) niove upwards (down- 
wards). On the other hand, the Ga atoms (both on surface and subsurface) are shifted 
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towards the bulk. 

In Intca, Mn relaxes inward by ~ 0.32 A; the Ga^nr/ and Gasubsurf atoms are displaced 
downwards while the Assurf i-Assubsurf) atom moves upwards (downwards). In summary, 
both in Intca and Int^^, cations (surface and subsurface) close to Mn move downwards, while 
anions upwards or downwards according whether they are on surface or subsurface. The net 
result is a local reduction of the surface buckling with respect to the clean unperturbated 
surface, more than 30 % and 40 % for IntAs and Intca respectively, with a net local buckling 
of about 0.46 A for Int^is and 0.40 A for Intca- As far as the interatomic distances between 
Mn and the nearest atoms are concerned (Tab.HJ, they are in general longer than the ideal 
bulk value by ~ 2-3 %; the distances between Mn and more distant atoms are shorter than 
the bulk cases, except for G^subsurf in IntAs? as it can be seen in Tab. HJ 

2. Electronic properties 

In fig.m we show the Density of States projected onto surface layer (PDOS); the continu- 
ous lines refer to Intyi^ or Intca while the dashed lines refer to the clean GaAs (110) surface. 
DOS for IntAs(IntGa) are shown to the left (right) side; the Fermi level (Ej) is set to zero 
eV. The d Mn projected DOS is also shown (grey area). The positive and negative DOS 
correspond to spin- up and spin-down components. First of all, in both Intyi^ and Intca, the 
DOS curves for Intyi^ and Int^a are very close to those corresponding to the clean surface 
case, but they differ in the energy region around Ef. An energy gap around is present 
in both majority and minority DOS. In Intyis the majority and minority spin gaps overlap 
and almost coincide, maintaining the surface semiconducting with a gap of about ~ 0.2 eV. 
In Intca, instead, majority and minority spin gaps are quite different: ~ 0.3 eV for the 
majority component and ~ 0.1 eV for the minority component. The perturbation is weak 
on the valence band and stronger on the conduction band. The main difference between 
Int^s and Intca DOS curve concerns a peak in the minority component in Intca around the 
Fermi energy (in IntAs it is shifted by 0.3-0.4 eV below the Fermi energy) which reduces the 
gap in Intca- 

In both systems, the Fermi level lies in the lower tail of the conduction band thus indi- 
cating that interstitial Mn impurity behaves as a donor, like in the bulk case.— 

At variance with the bulk case, where the calculated DOS for the two tetrahedral in- 
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terstitial positions are almost the same,— ^ thus indicating a week influence of the nearest 
neighbors on the interstitial Mn in the two configurations, the difference between surface 
Intca and IntAs cases is more sizeable, indicating a stronger effect of the local environment. 

The PDOS almost recover the bulk features already in the second layer (not shown in 
fig.01). Therefore, the introduction of Mn results in a perturbation of the electronic properties 
mostly localized on the first layer and strongly depending on the local environment. 

As far as the d states are concerned, we observe that their contribution to the occu- 
pied majority spin component is by far larger than their contribution to the minority spin. 
However, their overall weight in the GaMnAs system is negligible and the valence band is 
in practise almost non spin-polarized (as observed above). In both cases, the Mn spin- up 
d states are occupied and quite similar in shape while the spin-down d states are almost 
unoccupied and they have a different shape, especially around the Fermi level. 

In conclusion, the two Mn local environment give rise to a quite different surface electronic 
structure, with the differences mainly localized around the Fermi level. 

3. Magnetic properties 

In the following, we analyze the magnetic properties. The total and absolute magnetiza- 
tion in the supercell are different in the two configurations. They are equal to 4.23 and 4.84 
fiB in IntAs and to 3.41 and 4.71 /i^ in Intca- The difference between total and absolute 
magnetization corresponds to the presence of region of negative spin-density in the unit cell; 
this difference is higher in Int^a than in Int^s, suggesting higher (absolute) values and/or 
more extended region of negative spin-density in the former than in the latter. It also justi- 
fies the smaller total magnetization of Intca with respect to Inty^^- This is a clear evidence 
that the induced magnetization is strongly infiuenced by the local Mn environment. 

Interesting information can be gained by looking at the individual atomic magnetic mo- 
ments obtained as the difference between the calculated majority and minority Lowedin 
charges.— The results have been reported elsewhere.— The highest value of Mn spin- 
polarization is found in Int^s (3.96 /i^) while it is slightly lower in Int^a (3.67 /i^)- The 
Mn magnetic in ItAas is almost integer in agrement with the existence of a clear gap in the 
Mn-projected DOS and the unoccupied states just cutting the Fermi energy. It is worth 
noting that our calculated Mn magnetic moments are larger than those corresponding to 
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the Interstitial Mn in the bulk and they are rather close to the value indicated for ferro- 
magnetically coupled substitutional Mn impurities on the Ga sublattice in bulk GaAs. In 
fact, ab-initio calculations^L^L^ report a Mn magnetic moment for bulk Int^s equal to 2.70 
fiB- A recent experimental work^i show that Mn impurities on GaAs(llO) surfaces have 
magnetic moments significantly larger compared to the bulk case. The experimental and 
theoretical results would suggest in general an enhancement of the Mn magnetic moments 
due to surface effects. Our calculations, compared with previous bulk DFT studies , 
support this indications. 

For Int^s, the Assurf and Assubsurf atoms have a ferromagnetic coupling to Mn, with 
a small magnetic moment equal to 0.05 /i^. The induced polarization in more distant As 
atoms is totally negligible. The Ga^ur/ atoms couple antiferromagnetically with Mn with an 
induced polarization on it equal to -0.14 fiB- Other atomic moments are negligible. 

As far as the Intca configuration is concerned, a negative magnetic moment is induced 
on G^surf (-0.17 /is) while the Ga.subsurf atoms have a negligible polarization. The Assurf 
shows only a negligible polarization, while it is positive and equal to 0.05 /i^ for Assubsurf- 

Our results for the magnetic properties can be summarized as follows: in both cases, the 
cations couple antiferromagnetically to Mn spin moment while anions couple ferromagneti- 
cally. Furthermore, only surface cations are spin polarized, while both surface and subsurface 
anions do polarize. 

IV. STM IMAGES 

In figiniwe show the schematic front and side views of the relaxed underlying structure 
lattice and the XSTM images, for empty states at a reference positive bias voltages (+2.0 
V). In Intca, the two NN surface Ga atoms of Mn appear very bright with features extending 
towards the Mn and the atoms in the neighbourhood also look brighter than normal. For 
Int^s, a very bright elongated spot in the center of the surface unit cell delimited by As is 
visible. We would like to point out that the simulated XSTM images have clearly different 
shape for the two geometric configurations, so the two different local coordination should 
be distinguished by STM analysis. Further, the simulated STM images for ItAas case well 
compare with experimental XSTM images of the 5-doped samples.— 
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V. CONCLUSION 



In summary, we have used first-principles simulations to characterize Mn interstitial im- 
purity on the GaAs (110) surface. Strong local distortion on the (110) GaAs surface are 
introduced by Mn, especially when it is surrounded by Ga atoms. In both case, Mn polarizes 
the NN and NNN atoms, giving rise to a ferromagnetic Mn-As and to an antiferromagnetic 
Mn-Ga configuration. The simulated STM images show very different shape of the imaged 
Mn atom, suggesting that two configuration can be clearly differentiated by STM analysis. 
Finally, recent experimental STM images are qualitatively similar to our simulated one for 
Int^s configuration, suggesting the possible identification of Mn interstitials surrounded by 
As atoms in the experimental samples.— 
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TABLE I: Vertical atomic displacements with respect to ideal zinc blende bulk positions (first row) 
and nearest-neighbor surface and subsurface relaxed interatomic distances (second row) for Int^s 
(upper part) and Intca (lower part); +/— refer to an downward/upward relaxation; the numbers 
in round brackets refer to unrelaxed interatomic distances. Units are in A. 



Nearest Neighbor bond-lengths (A) 


IntAs 


ASsur/ 


■^^subsurf 


Ga^ur/ 


G&subsurf 


+0.15 


-0.19 


-0.06 


-0.06 


2.52(2.40) 


2.44(2.40) 


2.49(2.78) 


2.90(2.78) 






^^subsurf 




■^^subsurf 


-0.22 


-0.24 


+0.06 


-0.10 


2.48(2.40) 


2.56(2.40) 


2.68(2.78) 


2.63(2.78) 
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FIG. 1: Schematic side and top view of the clean GaAs(llO) surface. Only the three topmost layers 
(1** layer is the surface layer) are shown in the Figure. In this and other figures, black spheres are 
cations (Ga atoms), grey spheres are anions (As atoms). 
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FIG. 2: Conventional bulk unit cells representing Mn atom in tetrahedral-interstitial configurations, 
surrounded by As atoms (grey spheres) as nearest neighbors (top part, to the left) and by Ga atoms 
(black spheres) as nearest neighbors (top part, to the right). Bottom part: hexagonal interstitial 
position with Mn surrounded by 3 As and 3 Ga as nearest neighbors. 
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FIG. 3: Schematic side and top view of the relaxed Int^s and Intca configurations. Mn is explicitly 
shown. 
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FIG. 4: Density of States (DOS) projected into surface (continuous line) layer for Int^s (to the left) 
and Intca (to the right). Dashed line corresponds to the DOS for the clean surface. Mn projected 
DOS is also shown (grey filled area). The Fermi level is set to zero eV. 
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FIG. 5: Simulated STM images of isolated Mn interstitial in GaAs(llO) surface, with As NNs (to 
the left) and Ga NNs (to the right). Top panels: ball-and-stick model of the relaxed surface, top 
and side view (Ga: black spheres, As: grey spheres). Bottom panel: simulated STM images for 
positive bias voltage. The intersection of the dotted lines locates the position of Mn (projected on 
the (110) plane). 
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